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Abstract. We present a multi-wavelength study of the star forming region ISOSS J 20298+3559, which was identified by a 
cross-correlation of cold compact sources from the 170 fim ISOPHOT Serendipity Survey (ISOSS) database coinciding with 
objects detected by the MSX, 2MASS and IRAS infrared surveys. ISOSS J 20298+3559 is associated with a massive dark cloud 
complex (M ~ 760 Mq) and located in the Cygnus X giant molecular cloud. We derive a distance of 1800 pc on the basis of 
optical extinction data. The low average dust temperature (T ~ 16 K) and large mass (M ~ 120 Mq) of the dense inner part of 
the cloud, which has not been dispersed, indicates a recent begin of star formation. The youth of the region is supported by the 
early evolutionary stage of several pre- and protostellar objects discovered across the region: I) Two candidate Class objects 
with masses of 8 and 3.5 Mq, II) a gravitationally bound, cold (T ~ 12 K) and dense (n(H2) ~ 2 ■ 10 cm -3 ) cloud core with 
a mass of 50 Mq and III) a Herbig B2 star with a mass of 6.5 Mq and a bolometric luminosity of 2200 Lg, showing evidence 
for ongoing accretion and a stellar age of less than 40000 years. The dereddened SED of the Herbig star is well reproduced by 
an accretion disc + star model. The externally heated cold cloud core is a good candidate for a massive pre-protostellar object. 
The star formation efficiency in the central cloud region is about 14 %. 

Key words, stars: formation, accretion, Herbig Be, ISM: clouds, dust, ISM:individual objects: ISOSS J 20298+3559 



1. Introduction 

It is a challenge to identify massive young stellar objects 
during their early evolution. The youngest protostars form 
deeply embedded in their cold (T ~ 10-20 K) parental clouds 
(Pudritz 2002). The association with dense ambient material 
makes such objects best detectable as cold condensations at far- 
infrared and (sub)millimeter wavelengths. The short evolution- 
ary timescales (Palla & Stahler 1993) and low spatial density 
of massive objects require large scale surveys for their identi- 
fication. Many of the known intermediate- and high-mass pro- 
tostellar candidates have therefore been discovered by follow- 
up studies towards IRAS sources (eg. Shepherd et al. 2000, 
Cesaroni et al. 1997, Molinari et al. 1998, Beuther et. al 2002), 
which were selected on the basis of color and flux density cri- 
teria (e.g. by Wood & Churchwell 1989, Palla et al. 1991) 

The earliest stages of massive star formation are character- 
ized by the initial conditions of their parental cloud cores with 
spectral energy distributions peaking beyond 100 /im (Evans et 
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al. 2002). In order to unveil such young objects we are using 
the ISOPHOT (Lemke et al. 1996) 170 ^m Serendipity Survey 
(ISOSS) (Bogun et al. 1996), which is the largest high spatial 
resolution survey performed beyond the IRAS 100 /im band. 
We selected bright and compact sources detected by ISOSS and 
IRAS with a flux ratio [F170/im/F100/im] > 2, implying a low 
dust temperature T < 18 K and a large mass of the cold ISM in 
these objects. Since the clustered mode of massive star forma- 
tion commonly involves young stellar objects of different evo- 
lutionary stages, we require the presence of embedded sources 
with thermal infrared excess as indicated by the 2MASS (Cutri 
et al. 2000) and MSX (Price et al. 2001) infrared surveys. The 
latter criterium also avoids confusion with cold interstellar cir- 
rus. Here, we present the results of follow-up observations of 
the cold star forming region ISOSS J 20298+3559 and show 
evidence for its early evolutionary stage. 



2. Observations 

2. 1. ISO far-infrared measurements 

170 urn data covering the region around ISOSS J 20298+3559 
have been extracted from the ISOPHOT Serendipity Survey. 
The scans performed with the C200 2x2 pixel detector array of 
stressed Ge:Ga with a pixel size of 89.4 arcsec provide a spatial 
resolution of 1.8 arcmin FWHM. A flatfield correction for the 
individual scans was derived from redundant measurements at 
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Fig. 1. a) 170 fim continuum emission centered around ISOSS J 20298+3559, colors are from 180 to 400 MJy/sr. Overlaid are 
IRAS 100 /mi contours from the HIRES processing (contour level 240 MJy/sr). The map gives an impression of the typical sky 
coverage of the ISOPHOT Serendipity Survey towards the northern galactic plane, b) 170 /im surface brightness along an ISOSS 
scan (TDT16380300) across the region (arrow in Fig. [j]a), reconstructed from measurements with 2 pixels of the ISOPHOT C200 
camera. The brightness profile has been decomposed into a contribution by a compact source (FIR1 ) and an extended component. 



scan crossings. The calibration accuracy is estimated to be ± 
30 % (Stickel et al. 2000) 

2.2. SCUBA & MAMBO (sub)millimeter mapping 

Submillimeter jiggle maps at 450 /im and 850 /im were ob- 
tained with the SCUBA bolometer array (Holland et al. 1999) 
at the JCMT from Mauna Kea (Hawaii) on July 28, 2001 under 
excellent and stable sky conditions (rgso^m = 0.185 ± 0.01). 
The observing time was 30 min. The atmospheric transmission 
was determined from sky dips and water radiometer data (at 
JCMT and CSO). Mars, Uranus and CRL 618 served as calibra- 
tors. The data were reduced using the SCUBA User Reduction 
Facility (SURF) including identification of noisy bolometer 
pixels and removal of sky noise. The photometric accuracy 
derived from the calibration observations is 25 % at 450 /im 
and 20 % at 850 /im. The measured HPBW is 7.9 arcsec at 
450 /im and 14.9 arcsec at 850 /im (derived from observations 
of CRL618). In order to remove the error beam of the tele- 
scope, which contributes significantly to the 450 /im data, the 
SCUBA maps were deconvolved with a beam map of Uranus 
using the SCLEAN algorithm (Keel 1991) and finally restored 
to the initial spatial resolution. We know that the maps may be 
distorted by chopping into ambient cloud emission. This will 
affect the 450 /im flux more than the 850 /im flux, because the 
error lobe is larger at 450 /im. The 450 /im map was therefore 
finally corrected by adding a constant level of 5 % of the peak 
value to the map, before convolving it to the final resolution. 

1.3mm continuum observations were carried out on June 
9 and June 11, 2001 with the IRAM 30m telescope on Pico 
Veleta (Spain) using the 37-channel bolometer array MAMBO 
(Kreysa et al. 1998). Our observing strategy adopted On-Off 
measurements at 8 different raster positions in order to sam- 
ple most of the object plane. The mode was chosen due to a 



rather high sky opacity of ri.2mm = 0.62 and T\.2mm = 0.55 
and unstable atmosphere on the two days. Calibration measure- 
ments were obtained with the nearby source K30A. Data were 
reduced with the NIC package including sky noise reduction. 
The photometric accuracy is 30 %. 

2.3. UKIRT/MAX mid-infrared imaging 

Diffraction limited mid-infrared images were acquired with the 
MAX camera mounted at the 3.8m United Kingdom Infrared 
Telescope UKIRT from Mauna Kea (Hawaii) on December 20, 
2000. The MAX camera (Robberto et al. 1998) consists of a 
Rockwell 128x128 Si:As BIB array, providing a field of view 
of 35 x 35 arcsec 2 . Spectrophotometric data were obtained 
through three narrow band filters (A e // ~ 8.7, 9.7 11.6 /im, 
A A ~ 1.0 /im) and in the M-, N- and Q-bands. The obser- 
vations were made using standard beam switching techniques 
with a chop frequency of 2.5 Hz and a throw of 25 arcsec (PA 
90°). a Tau and a Boo served as standard stars, also for the nar- 
row band filters where we derived synthetic photometry from 
data obtained by Cohen (1995). The images were finally decon- 
volved by the algorithm from Richardson & Lucy with standard 
stars as PSF estimators and have a resolution with a FWHM of 
~ 0.6 arcsec. The photometric accuracy is 20 %. 

2.4. Optical spectroscopy and photometry 

Optical long-slit spectroscopy was obtained using the TWIN 
double spectrograph at the 3.5m telescope at Calar Alto (Spain) 
on December 12, 2000 with a spectral dispersion of 0.55 
A/pixel for the blue (3400-5500 A) and 1.1 A/pixel in the 
red (6300-7100 A) image channel. A photometric accuracy of 
10 % was obtained. The spectra were extracted according to the 
method of Home (1986). 
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Fig. 2. a) 850 pan continuum map of the compact ISOPHOT Serendipity Survey source FIR1, overlaid on a near-infrared JHK5- 
composite constructed from 2MASS data. Three compact dust condensations (SMM1, SMM2 and SMM3) are detected, which 
are located in a diffuse extended emission. Mid-infrared sources detected by the MSX-satellite are marked with boxes. The 
submillimeter knot SMM2 is associated with a small cluster of embedded NIR sources (IRS2..5) as detected by the 2MASS and 
MSX surveys. IRS 1 was identified as a very young Herbig B2 star by our follow-up spectroscopy. The two compact submillimeter 
sources SMM1 & SMM3 without any infrared counterparts are candidate Class objects. The position of the very cold cloud 
core is indicated by an ellipse. Contour levels are starting at 67 and increasing by 33 mJy/beam. The size of the SCUBA beam 
is indicated in the lower right, the dashed circle corresponds to the ISOPHOT beam, b) The dust color temperature distribution 
across FIR1 shows the presence of a very cold (T^ ~ 11 K) core at the center of the cloud. The temperature profiles towards 
north-east and south-west indicate an external heating by the infrared sources IRS 1-5. The temperature is calculated from the 
submillimeter spectral index between 450 /im and 850 /im, assuming a dust emissivity (3=2. 



CCD images in the Johnson BVRI bands were obtained 
with the MPIA 0.7m telescope on mount Konigstuhl in 
Heidelberg (Germany), equipped with a Tektronix 1024x1024 
camera (scale 0.89 arcsec/pixel). Data reduction of both spec- 
troscopy and direct imaging was carried out using standard 
IRAF routines. I-band photometry of about 3000 stars has been 
derived using DAOPHOT (Stetson et. al 1987). The stellar pho- 
tometry is complete to a limiting magnitude I = 18.8 mag with 
a photometric uncertainty a < 0.2 mag. 

2.5. Effelsberg NH 3 measurements 

We carried out NH3 measurements with the MPIfR 100m tele- 
scope at Effelsberg in April 2002. The HPBW at the observ- 
ing frequency of 23.7 GHz is 40". The facility 1.3 cm receiver 
was used with a typical system temperature (antenna tempera- 
ture units) on the sky of 170 K in total power mode. The back- 
end was the facility 8096 channel autocorrelator split into 4 
bands in order to observe simultaneously both polarizations at 
the NH 3 (1,1) and (2,2) rest frequencies of 23694.495 MHz and 
23722.633 MHz, respectively. The resulting spectral resolution 
was 0.06 km s _1 . Pointing and calibration were checked by 
continuum scans across NGC7027. We estimate a pointing ac- 
curacy of about 5". We calibrated our data assuming a main 
beam brightness temperature for NGC7027 of 8.2 K corre- 



sponding to 5.86 Jy (Baars et al. 1977). Our NH3 spectra were 
reduced using the GILDAS software, and analyzed following 
procedures described by Harju et al. (1993). 

3. Results 

3. 1. Morphology of the cold ISOSS source 

Fig. [I] presents our far-infrared measurements of ISOSS J 
20298+3559 obtained by the ISOPHOT Serendipity Survey. 
The source was detected by two individual scans of the sur- 
vey. A 170 /im continuum map constructed from all scans 
in the region (Fig. la) shows the compact object on a bright 
galactic background. Overlaid on the map are IRAS 100 /im 
contours, indicating the previous detection of the associated 
source IRAS20278+3549. Due to the superior spatial resolu- 
tion of ISOPHOT we could disentangle two emission com- 
ponents from a central scan across the object (Fig. [l|b). The 
170 /im emission is well described by an unresolved source 
(FIR1) with a flux density of 8170^771=1 10 Jy and an extended 
emission component (size 5.9 arcmin FWHM, peak surface 
brightness Iiro^m = 90 MJy/sr). 

High resolution sub-millimeter continuum maps at 450 /im 
and 850 fim with SCUBA and at 1.3 mm with MAMBO have 
been obtained towards the unresolved far-infrared source FIR1 . 
The 850 /im image superimposed on a near-infrared (JHK5) 
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Table 1. Photometry of (sub)millimeter sources associated with ISOSS J 20298+3559-FIR1 



Source 


Position [J2000] 






Total Flux 




Size 




Peak Flux 






a 


8 




Jy 




pc 




T 1 1 

mJy beam 










450 fim 


850 


1.2 mm 




450 /im 


850 


1.2 mm 


SMM1 


20:29:49.9 


+35:59:15 


1.26 


0.233 


0.075 


0.1 x 0.2 


690 


251 


80 


SMM2 


20:29:45.6 


+35:58:46 


0.31 


0.091 




<0.07 


610 


206 




SMM3 


20:29:48.2 


+35:59:24 


1.05 


0.155 




0.15 x 0.2 


290 


117 




cold core 


20:29:48.4 


+35:58:57 


6.02 


0.902 


0.33 


0.55 x 0.29 


410 


203 


60 


extended 






18.86 


2.47 


0.7 


0.9 x 0.7 


368 


173 





color composite from the 2 MASS survey is presented in Fig. 
||. FIR1 indicated by the dashed circle is clearly resolved by the 
JCMT in the submillimeter. The map indicates the presence of 
a dense core (outlined by an ellipse) surrounded by extended 
emission, which covers an area of about 100 x 70 arcsec 2 and 
accounts for the bulk of the flux. This extended emission is 
located northeast and southwest of the core and shows three 
compact dust condensations. We derived positions and fluxes 
of the compact sources by a two-dimensional gaussian decon- 
volution using a local background estimation. The source sizes 
were determined from the 450 /im data, which have the highest 
available resolution. The 1.2 mm fluxes of the fainter sources 
SMM2 and SMM3 could not be determined due to background 
confusion at that wavelength. The results are given in Tab. [j]. 
The integrated (sub)millimeter flux density of FIR1 within a 
circular 2.4 arcmin aperture corresponding to the SCUBA field 
of view and covering most of the ISOPHOT beam are S45o A tm 
=27.5 Jy, S 8 5o M m = 3.9 Jy and S^oo^m = 1.1 Jy. 

3.2. Distance and structure of the cloud complex 

A well known drawback for the use of kinematical distances 
towards Cygnus is the near-zero radial velocity gradient up to 
a distance of about 2 kpc from the sun. Therefore we have 
created an extinction map to determine a reliable distance to 
ISOSS J 20298+3559 based on a Wolf diagram and to study the 
large-scale distribution of dust associated with the far-infrared 
source. We have obtained deep I-band images of the region us- 
ing the MPIA 0.7m telescope. 

The extinction map shown in Fig. [| was constructed 
according to the method of Bok (1956) and subsequently 
smoothed to the ISOPHOT resolution. An adaptive grid (re- 
seau size ~ 1 arcmin) was used for optimization of the angular 
resolution with respect to the stellar density (Cambresy et al. 
1997). The observed cloud extinction A obs with respect to the 
reference positions was calculated from the number density N b 
of stars seen towards the cloud and their unobscured reference 
density N b ef : A obs = -1/a \og(H b ref /N b ), where a denotes 
the slope of the I-band luminosity function. We use Ay = 2.07 
Ai (Rieke & Lebofsky 1985) to relate visual extinction to our 
I-band observations. The reference fields are marked with 'R' 
in Fig. |]and are located in direction of the lowest dust column 
density as indicated by the 170 /im continuum. 

The extinction A obs derived from the simple Bok for- 
mula underestimates the true cloud extinction A true of a 
distant cloud due to the contamination by foreground stars. 




20:30:00 45 30 

Right ascension (J2000.0) 



Fig. 3. Extinction map towards the star forming region ISOSS 
J 20298+3559 derived from I-band star counts and corrected 
for foreground stars . Four individual dark clouds (C 1 . . C4) have 
been detected. Overlaid are contours of visual extinction Ay 
(contour levels at 1.5, 2.5, 3.5 and 4.5 mag). The rectangles 
marked with 'R' show the reference fields used to set the zero 
point of extinction. In order to obtain the total line-of-sight ex- 
tinction, a foreground extinction of Ay = 2.7 mag has to be 
added to the data. 

Foreground stars lower the density contrast towards the cloud. 
Since they are not affected by the cloud extinction, their num- 
ber density remains constant: = Nl In order to account 

re. J 

for the stellar foreground and background population we apply 
the general form of the Bok formula (e.g. Cambresy et al. 2002) 
providing the true extinction A 1 /^ depending on the fraction 
of background stars. 

Fig. ^presents the Wolf diagram for the reference field (up- 
per curve) and towards the dense region {A obs > 1 mag) of the 
cloud complex (lower curve). The logarithm of the cumulative 
distribution of apparent magnitudes is shown, error bars are ± 1 
er and are calculated according to Poisson statistics. The num- 
ber density in the reference field agrees well with the predic- 
tion from the Galactic stellar distribution model from Robin & 
Creeze et al. (1986), assuming an average interstellar extinc- 
tion coefficient of Ay = 1.5 mag kpc -1 in the distance range 
< d < 4 kpc and Ay = 0.7 mag kpc -1 beyond 4 kpc. The 
higher extinction below 4 kpc is attributed to dust in several 
giant molecular cloud complexes of the local spiral arm and 
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Fig. 4. Wolf-diagram of cumulative I-band star counts towards 
the dark cloud (lower points) and the reference fields (upper 
points). The error bars correspond to ±lcr of a Poisson statis- 
tics. The data are in good agreement with our predictions based 
on a Galactic stellar distribution model, indicated by the over- 
laid curves. The results of the Monte-Carlo simulations are 
consistent with a distance of 1.8 ± 0.3 kpc to the dark cloud 
complex 



extending out to the inter-arm gap between local and Perseus 
spiral arms at a distance of 4 kpc (Ungerechts et al. 2000). The 
integrated (0 < d < 5 kpc) line-of-sight extinction of Ay = 
6.7 mag towards the reference fields agrees with the value of 
Ay = 6.9 ± 0.7 mag derived by Schlegel & Finkbeiner (1998) 
from the IRAS 100 /tm dust continuum. 

To interprete the shape of the Wolf diagram correctly, 
a model curve of an obscured region was determined from 
Monte-Carlo simulations using a Galactic stellar distribution 
model for the unobscured field star luminosity function. We as- 
sumed the extinction in the cloud to be gaussian distributed for 
resolutions smaller than tested by the average extinction AJ fcs 
of our integral aperture. The lower curves in Fig. ^ show the re- 
sults of our simulations for three best-fit cloud models, which 
are consistent with our observed average extinction (Aj bs = 1.1 
mag) and are located at distances of 1.5, 1.8 and 2.1 kpc, re- 
spectively. From the observations we derive a cloud distance d 
= 1.8 ± 0.3 kpc and an average extinction of the cloud of Ay 
= 3.1 ± 0.5 mag, where the errors mainly originate from the 
Poisson statistics of our star counts. Using the interstellar ex- 
tinction coefficient of Ay =1.5 mag kpc -1 determined from 
the reference field, the foreground extinction towards ISOSS J 
20298+3559 is Ay = 2.7 ± 0.3 mag. 

We compared our result with the empirical distance esti- 
mator d = 320 N 57 pc by Herbst & Saywer (1987), which 
relates the cloud distance to the number N of foreground stars 
counted within a circular aperture of 5 arcmin diameter on the 
POSS blue print. Using B-band star counts from the USNO- 
PMM catalogue (Monet 1996) we obtain a distance of d ~ 2.2 
kpc for the cloud complex. The formula is however only valid 
for opaque dark clouds (Ay > 7 mag) and the higher distance 



Table 2. Individual clouds towards ISOSS J 20298+3559. The 
columns are: (1) cloud, (2) position of peak extinction Ay, (3) 
mass Mm + M# 2 , (4) peak extinction Ay and (5) size. 



Source 


Position [J2000] 


Mass 


Peak Av 


Size 




Q, S 


M Q 


mag 


pc 2 


CI 


20:29:46.8 +35:59:30 


188 


4.7 


2.0 x 1.7 


C2 


20:29:35.1 +35:58:09 


162 


4.6 


1.9 x 1.7 


C3 


20:30:05.3 +36:00:07 


115 


2.6 


2.5 x 1.5 


C4 


20:29:31.9 +36:00:49 


295 


3.7 


2.3 x 2.2 



might originate in the partly translucent character of ISOSS J 
20298+3559. 

3.3. Cloud masses 

The extinction map clearly reveals the presence of a dusty 
cloud complex with 4 different dark clouds (C1..C4). CI coin- 
cides with the cold source FIR1 and accounts for the extended 
far-infrared emission shown in Fig. [fjb. Despite being covered 
by our 170 /im scan, the south-western cloud C2 does not show 
up with a strong 170 /im excess and we find an upper limit of 
Si70/jm < 15 Jy. The cloud complex also shhows a small ex- 
tinction ridge towards the east (C3) and extends to the north- 
west (C4). 

Cloud masses were derived from the extinction map, using 
the relation by Dickman (1978): 

M = (Add) 2 fJ--j^- y~] A v (i), (1) 

Ay '—' 

i 

where A57 is the angular size of the map, d is the distance to 
the cloud, /i is the mean molecular weight corrected for helium 
abundance, and i corresponds to a map pixel. With the dust- 
to-gas ratio proposed by Savage & Mathis (1979), Nn/Ay = 
1.87 x 10 21 cm" 2 mag" 1 (N H = N HI + 2N H2 ) we obtain a 
total mass of the cloud complex of 760 M Q . We list in Tab. 
H parameters of the four individual clouds. Due to systematic 
uncertainties in the extinction values and the conversion factor, 
we assume the mass estimates to be accurate within a factor of 
2. 

3.4. Dust properties of FIR1 

Fig. U displays the spectral energy distribution (SED) for the 
cold region FIR 1 . The spectral energy distribution longwards 
of 100 /im is dominated by large grains, which are at an equi- 
librium temperature within the surrounding radiation field. 

The flux density of the diffuse dust component in FIR1 (ex- 
cluding the knots SMM1, SMM2 and SMM3) is well character- 
ized by optically thin thermal radiation of two modified black- 
bodies Si, oc yPB(y, Td) with a dust emissivity (3=2 and dust 
temperatures Td = 13 ± 2 K and T^ = 19 ±2 K. As can be seen 
from Tab. [TJ the contribution of SMM1, SMM2 and SMM3 to 
the total (sub)millimeter flux is less than 15 %. In the large 
beam measurements of the IRAS/ISO satellites we cannot sep- 
arate these sources. If they contributed significantly more FIR 
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Table 3. Results of the ammonia observations of the cold core region. The columns are (1) line, (2) main beam brightness 
temperature of the maingroup, (3) Vlsr, (4) FWHM from fit, (5) main group line area, (6) total line area of satellite groups and 



Line 


T M B(main) 


^LSR 


FWHM 


areamain 


areasatel 


r 




K 


kms -1 


kms~ 


Kkms -1 


Kkms- 1 




(2,2) 


0.139 (0.074) 


0.66(0.13) 


0.883(0.408) 


0.146 




0.78(0.10) 


(1,1) 


0.329 (0.076) 


0.58(0.02) 


0.955(0.065) 


0.34 


0.19 


0.10(0.05) 




Wavelenglh [^m] 

Fig. 5. Spectral energy distribution of the cold ISOSS source 
FIR1. The total flux density is well characterized by the super- 
position of optically thin thermal radiation from two modified 
blackbodies (j3 = 2) with dust temperatures of T^ = 13 ± 2 K 
and Td = 19 ±2 K, respectively. 

flux than estimated from our submillimeter data, the dust tem- 
peratures of the diffuse dust component would be even lower 
and the resulting dust mass higher. 

Drawing the conclusion of a very cold dust component 
merely from large-beam integrated photometry is critical with- 
out any radiative transfer model. Our spatially resolved sub- 
millimeter data allow, however, to confirm the presence of a 
very cold dust component from the dust color temperature map 
of ISOSS J 20298+3559-FIR1 presented in Fig. |b. The map 
was obtained from the deconvolved 450 /im and 850 fim maps 
smoothed to 14.9 arcsec resolution. The color temperature T co / 
was derived from the spectral index between the two wave- 
lengths and is given by 



l)/vl +l3 (z hul/kTc ° l - 1). ( 2 ) 



with a dust emissivity (3 = 2. A decreasing color temperature 
and a clear 850 /im (and 1.2 mm) excess is observed towards 
the cloud center, exhibiting a minimum temperature of only 
T co/ = ll±2K. 

Optical depth effects in the sub-millimeter or a change in 
the dust emissivity (j3 < 2) could affect the spectral index. 
We conclude, however, that these effects are small and the dust 
temperature dominates the spectral index: The peak column 
density towards the cloud core derived from our 850 /im data is 
NCH2) = 2.2 • 10 22 cm -2 , in good agreement with our ammo- 




Velocily (km/s) 



Fig. 6. Effelsberg NH3 spectrum taken at the position of the 
cold core region (Table 1) in ISOSS J 20298+3559-FIR1 



nia observation (see next section). The corresponding optical 
depth is Tg5o = 2 • 10 -3 << 1. Assuming a roughly ellip- 
soidal core geometry of 0.55 x 0.29 pc 2 as indicated in Fig. 
we obtain an absolute molecular gas density of N(H2) = 2 • 10 
cm -3 . Concerning j3, the centrally decreasing temperature pro- 
file points to an external heating of the core region. Therefore 
the conditions in our core region basically resemble observa- 
tions of nearby very cold pre-stellar cores for which (3=2 was 
found (eg. Juvela et al. 2001, Ward-Thompson et. al. 2002). 

Masses for the emitting regions have been derived from 
the optically thin emission in the submillimeter. The dust mass 
is given by M d = F S50 d 2 / (K 850 B 850 (T d )), where we used a 
dust mass absorption coefficient of k 8 5q =2.1 cm 2 g _1 follow- 
ing Ossenkopf & Henning (1997) for dense protostellar cores. 
Following Sodroski et al. (1997) we assume a gas-to-dust mass 
ratio of 150 and derive a total gas mass M/f 2 =120 M© ± 24 
M Q for FIR1. The total mass contained in the cold core is Mjj 2 
= 50±10M o . 

We calculated the bolometric FIR luminosities re-radiated 
from large grains by integrating over the Planckian radiation 
and find for the region Lpjn = 250 L© and for the cold core 
Lfir =15 L0, respectively. 

3.5. Dense gas in the cold core region 

It is important to approve the physical conditions derived from 
our dust continuum observations by molecular line observa- 
tions, which reflect the physical conditions of the gas phase 
ISM. NH3 is a very useful indicator of molecular cloud tem- 
peratures and column densities (e.g., Walmsley & Ungerechts 
1983). 
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Fig. 7. a) 8.3 /im image obtained by the MSX satellite with contours of visual extinction overlaid (see Fig. [3[ contour levels 
Ay=1.5, 2.5, 3.5 and 4.5 mag). Colors are stretching from 5 to 40 mJy pixel -1 . Five sources are detected by MSX above the 5 
a level, the two fainter ones are marked by arrows. Note the extended source towards the south-east, which shows an external 
illumination of the outer rim of the molecular cloud. The illuminating star is IRS7. b) I-Band image of the ISOSS J 20298+3559 
region, obtained with the MPIA 0.7m telescope with contours of visual extinction overlaid. The optical identifications of mid- 
infrared sources seen by the MSX satellite are labeled. 



We detected the NH 3 (1,1) line with a S/N of 4 and the 
NH 3 (2,2) fine with a S/N of 2. Fig. | presents the (1,1) spec- 
trum. We derived the NH 3 (1,1) line optical depth r n = 0.10 ± 
0.05 and estimated the (1,1)(2,2) rotational temperature T 12 = 
19.2 ± 4.0 K and the gas kinetic temperature T kin = 21.3 ± 5.0 
K. The derived ammonia column density was estimated as 
N(NH 3 )=1.32±0.24 x 10 14 cm- 2 . From the AT(H 2 ) column 
density obtained by our submillimeter continuum observations 
and smoothed to 40 arcsec resolution we derive a molecular hy- 
drogen column density of N(H 2 ) = 1.3 ± 0.6 x 10 22 cm -2 for 
the region covered by the Effelsberg beam. This corresponds 
to a fractional abundance of x( NH 3 ) = N(NH 3 )/N(H 2 ) = 
1.0- 10" 8 which is in the range found by Molinari et al. (2000) 
in their study of ammonia clumps associated with young in- 
termediate and high mass star forming regions. Ammonia in 
ISOSS J 20298+3559 seems however to be under-abundant 
when comparing to the values given by Harju et al. (1993). 
They find a mean value of x(NHz) ~ 2.6 • 10~ 8 in a sam- 
ple of 22 ammonia clumps in Orion. Ammonia is considered 
as a molecule characterizing later stages of chemical evolution. 
Myers & Benson (1993) found NH3 to be more abundant in 
older cores, where stars have already formed. This may indi- 
cate that ISOSS J 20298+3559 is a rather young object. Tab. | 
summarizes our ammonia results. 

ISOSS J 20298+3559 is located towards the southern out- 
skirts of the Cygnus X Giant Molecular Cloud (Dame & 
Thaddeus 1985). The majority of active and luminous star 



forming regions in this giant molecular cloud is found at dis- 
tances between 1.2 and 2.5 kpc (Odenwald & Schwartz 1993), 
corresponding to radial velocities (LSR) between -10 km s _1 
and +10 km s — 1 , The ammonia line velocity Vlsr = 0.58 ± 
0.06 km s _1 is consistent with the kinematical association with 
the Cygnus X GMC. 

3.6. Young stellar objects in the cloud complex 

Mid-infrared observations are well suited to identify young 
stellar objects enshrouded by warm circumstellar dust. An 
MSX 8.3 /im image covering the region of our extinction map 
is presented in Fig. Four compact and one extended mid- 
infrared sources have been detected above the 5 a noise level (a 
~ 15 mJy/beam). The astrometric accuracy of the MSX image 
(~ 2 arcsec RMS) allows the cross-identification with objects 
from our I-band image shown in Fig. ^|b and the 2MASS point 
source catalogue. All four compact MSX sources coincide in 
position with very red (I-K5 = 5. .7 mag) near-infrared coun- 
terparts. The objects are marked in Figs. ^ and [7]b. Astrometry 
and photometry of the sources are given in Tab. where the 
8.3 /im fluxes were derived from the MSX map performing 
aperture photometry, yielding a photometric accuracy of 25 %. 

The intrinsic excess emission and extinction of the MSX 
near-infrared counterparts are characterized by the (J-H) vs. 
(H-K s ) color-color diagram presented in Fig. ||, where the col- 
ors of IRS1..8 are plotted with their photometric 1 a errors. 
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Table 4. Infrared sources associated with the cloud complex The columns are: (1) IR source label, (2) position (3-6) near-infrared 
brightness in IJHK S bands, (7) 8.3 /jm flux, (8) 2MASS PSC identification, (9) additional remarks. 



Source 


Position [12000] 


I 


I 


H 




8.3 /J,m 


2MASS 


remarks 




a, 5 




mag 




mly 


counterpart 




IRS1 


20:29:47.9 +35:59:26 


12.90 


10.98 


9.45 


8.15 


1740 


12029479+355926 


Herbig B2 star 


IRS2 


20:29:45.8 +35:58:47 


17.23 


14.70 


13.01 


11.95 


90 


12029458+355847 


brightest member of cluster IRS2..5 


IRS 3 


20:29:46.1 +35:58:46 


18.54 


16.1* 


14.3* 


13.1* 






IHK S fluxes from 2MASS map 


IRS4 


20:29:46.1 +36:58:52 


19.5 


16.2* 


13.85 


12.25 




12029461+365852 


I flux from 2MASS map 


IRS5 


20:29:46.4 +35:58:50 


18.44 


15.74 


13.85 


13.02 




12029464+355850 




IRS6 


20:29:36.0 +35:58:18 


15.77 


13.33 


11.83 


10.73 


75 


12029360+355818 




IRS7 


20:29:54.5 +35:57:38 


12.48 


11.38 


10.97 


10.74 


1870 


12029545+355738 


MIR nebula + illuminating star 


IRS 8 


20:29:47.9 +36:01:08 


13.09 


9.78 


7.94 


6.98 


225 


12029479+360108 


background star 



7.0 - 



i e' 

L-IRS5 X 
Xf- IRS8 X 

S\ ± s' 




Fig. 8. 2MASS JHK5 color-color diagram towards the cloud 
complex. The number density of field stars is indicated by con- 
tours. The solid curve in the lower left of the diagram is the 
locus of unreddened main sequence stars and giants, their red- 
dening band is indicated between the two dashed lines. Several 
of the MSX counterparts IRS1..8 show an infrared excess char- 
acteristic for young stellar objects, as indicated by the colors 
of Herbig Ae/Be and T Tauri stars. Also given is the redenning 
vector from Rieke & Lebofsky (1985). 



The solid curve corresponds to the loci of unreddened main 
sequence and giant branch stars (Bessel & Brett 1988) and is 
extended by the reddening band, which confines stars with nor- 
mal photospheres (dashed parallel lines). 

IRS1, IRS2, IRS3, IRS4, and IRS6 lie to the right side 
and show an intrinsic infrared excess due to warm circum- 
stellar dust. Their observed colors are consistent with red- 
dened (Ay ^ 15 mag) pre-main-sequence stars of low- to 
intermediate-mass as indicated in the diagram by the loci of 
unreddened Herbig Ae/Be (dashed parallelogram, Hillenbrand 
et al. 1992) and T Tauri stars (dashed line, Meyer et al. 1997). 
The derived extinction is smaller than the total line-of-sight ex- 
tinction derived from our 850 /im, yielding Ay=18 mag to- 
wards IRS1 and Ay =23 mag towards the small cluster IRS2, 
IRS3, IRS4 and IRS5. The sources are not affected by the to- 
tal dust column, and are therefore either foreground objects or 



embedded in the cloud complex. The positional coincidence of 
dust with IRS1..5 (Fig. ||) favors however the association with 
the dark cloud: IRS1 is located on a dust ridge in the north- 
west of the central region of the cloud CI, IRS2 coincides in 
position with the compact submillimeter source SMM2 and is 
surrounded by a compact cluster of a least 3 further red sources, 
IRS5, IRS7 and IRS8 show no infrared excess. IRS5 and 
IRS8 are reddened with about (Ay ~ 15 mag). The corre- 
sponding line-of-sight extinctions from our extinction maps are 
23 mag for IRS5 and 4.5 mag for IRS7 and IRS8. which are lo- 
cated at the border of the cloud. This suggests that also IRS5 
and IRS7 are associated with the cloud, while IRS8 is a back- 
ground giant. 

IRS7 is located approximately in the center of a spherically 
segmented mid-infrared nebula seen in Fig. Qa. The nebula is 
located in the south-east of the dark cloud CI and indicates an 
external illumination of the outer rim of the molecular cloud by 
IRS7. The nebula and IRS7 coincide in position with the source 
IRAS20279+3547. The flux ratio F12/F25 > 1 .4 and the strong 
emission in the MSX 8.3 /im band, which contains a strong 
PAH feature, indicate a reflection nebula around a B type star. 
In order to clarify whether IRS7 is the illuminating source we 
have obtained B VRI photometry of the source. Comparing with 
the colors of main sequence stars we obtain a best fit for Ay = 
4.2 ± 0.3 mag and a spectral type B4..B6, in agreement with 
the total line-of-sight extinction Ay = 4.2 ± 0.5 mag towards 
IRS7 from our extinction map. 

3. 7. The embedded Herbig B2 star IRS 1 

In order to establish the pre-main-sequence nature of the bright 
source IRS1 we have obtained optical long-slit spectra pre- 
sented in Fig. ^[ The spectra display strong emission lines of 
HI, Fell, NI, Fel which are all characteristic for luminous 
young stellar objects (Cohen & Kuhui 1979). Parameters of im- 
portant spectral features are given in Tab. || 

We have several strong indications that IRS 1 is a rare ex- 
ample of an early B type pre-main-sequence object: He I ab- 
sorption lines of presumably photospheric origin have been de- 
tected at 6678 A and 7065 A. Comparing their strength with 
the theoretical values for stellar atmospheres given by Auer & 
Mihalas (1973) we derive a spectral type of B2 (T e// = 20000 
± 2500 K, log g = 4.0). This is consistent with the presence of 
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Fig. 9. Optical spectrum of the YSO ISOSS J 20298+3559-IRS1 taken with the TWIN spectrograph at the 3.5m telescope at 
Calar Alto. A spectral type B2 was inferred from photospheric He I absorption features. The spectrum is dominated by intense 
emission lines. They originate in a dense emission region, indicated by Fe II emission. 



Table 5. Observed strong features in the optical spectrum of 
ISOSS J 20298+3559-IRS1. The columns are: (1) transition, 
(2) observed peak wavelength, (3) equivalent width, (4) line 
flux. 



Table 6. Optical and mid-infrared photometry of ISOSS J 
20298+3549-IRS1 given in magnitudes. 



Line 


Obs. Wavelength 


EW 


Flux 




[A] 


[A] 


[10 -16 erg cm -2 s" 1 ] 


Ha 


6562.8 


-200 


1950 


H/3 


4861.4 


-19.9 


99.8 


Hy 


4340.4 


-5.2 


15.2 


Hel 


6678.15 


0.55 


-13.8 


Hel 


7065.18 


0.29 


-8.5 


Fell 


4923.06 


-2.8 


19.2 


Fell 


5017.36 


-3.5 


19.2 


Fell 


5316.05 


-4.3 


48.3 


Fell 


5167.75 


-3.4 


33.4 


Fel 


6358.56 


-0.2 


4.8 


Mgl(+Fel) 


6317.84 


-1.3 


26.8 


DIB 


6613.81 


0.3 


-6.5 


NI 


6442.02 


-0.4 


9.4 


NI 


6491.36 


-0.8 


16.9 


NI 


6505.58 


-0.8 


18.1 



NI emission lines which are confined to hotter Be type stars due 
to their high-ionization energy (14.5 eV) and the large line ex- 
citation 10 eV) (Andrillat et al. 1988). The relative strength 
of Fell with respect to neutral Fel also requires the presence of 
a more ionized line-emitting region than associated with clas- 
sical T-Tauri stars (Hamann & Persson 1992), where Fel lines 
often rival the strength of Fell. In order to derive the rotational 
velocity of IRS 1 from the He I lines we have used the rotational 
broadening function given by Gray (1976). We obtain v sin(i) 
= 158 ± 25 km/s, in agreement with a pre-main-sequence star 
of intermediate mass (Finkenzeller 1985). A small nebulosity 
(size 15 x 8 arcsec 2 , PA = 40°) around IRS 1 which is visible on 
POSS-B and -R plates was covered by our long slit spectra and 
identified as a reflection nebula. Its presence and the associa- 
tion with the dark cloud CI make the emission line object IRS 1 



Band 


B 


V 


R 


M 


P8 


S9 


N 


Pll 


Q 


Flux 

a 


17.3 
0.1 


15.7 
0.1 


14.5 
0.1 


5.7 
0.2 


3.9 
0.2 


4.1 
0.2 


3.7 
0.2 


3.1 
0.2 


2.6 
0.2 



by definition a Herbig Be star (Herbig 1960). According to the 
most comprehensive catalogue of Herbig Ae/Be stars compiled 
by The et al. (1994), only 20 objects known to date have a spec- 
tral type B2 or earlier. 

IRS1 (=MSX5C.G075.5314-01.8159) was detected by 
MSX at 8.3, 12.1 and 14.7 ^m and is included in the MSX 
point source catalogue (Egan et al. 1999). The source is un- 
resolved by the small MSX telescope (beam ~ 12 arcsec). In 
order to explore the nature of the thermal excess emission of 
warm dust around IRS1 we have obtained diffraction-limited 
sub-arcsecond images with the thermal infrared camera MAX 
at UKIRT. The observations cover the M, N and Q bands as 
well as narrow-bands centered on the 9.7 /m silicate feature 
and two strong PAH bands accessible from ground at 8.7 and 
11.6 fim. 

The high-resolution data confirm the presence of a com- 
pact mid-infrared source coinciding with IRS 1 . An unresolved 
point source (size < 0.5 arcsec FWHM) accounts for the bulk 
( ~ 80 %) of emission. The image obtained in the 11.6 /im 
PAH band shown in Fig. |l0| shows the presence of a faint ex- 
tended emission region towards the north of the stellar source. 
Optical broadband photometry of the source in the Johnson 
B VRI bands have been performed in order to cover the spectral 
energy distribution of the stellar photosphere of IRS 1 . Tab. ^ 
presents our optical and mid-infrared photometry of IRS 1 . 
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Table 7. Physical properties of the cold core in ISOSS J 
20298+3559. 



-5 5 

R.A. offset [a re sec] 

Fig. 10. Diffraction-limited 11.6 /im image of the Herbig B2 
star IRS1, obtained with the thermal infrared camera MAX at 
UKIRT using a narrow PAH-filter. The bulk of the emission 
originates from an unresolved compact source (size < 0.5 arc- 
sec FWHM). A faint region of extended PAH emission is seen 
towards the north. An image of the stellar point source HR 8316 
is inserted. 



4. Discussion 

4. 1. Evidence for a pre-protostellar core in the 
center of the dark cloud 

Our submillimeter observations towards the central region of 
the dark cloud complex have revealed the presence of a cen- 
trally condensed, dense (n(H2) = 2 • 10 5 cm -3 ), very cold 
(Trf ~ 12 K) and compact (R ~ 0.2 pc) core with a decreas- 
ing temperature gradient towards the center. Similar conditions 
have been recently observed in FIR studies of nearby pre- 
protostellar cores (Juvela et al. 2001, Ward-Thompson et al. 
2002) and are considered as the initial conditions which pertain 
in clouds from which stars form. Most of the known prestellar 
cores are however of low-mass. Is it possible that FIR1 contains 
a rare example of a massive pre-protostellar core. Even consid- 
ering that probably more than 50 % of the mass are dispersed 
during the early phases of star formation, the total mass of the 
core M = 50 M Q is likely to be massive enough to form at least 
one late O- or early B-type star. 

By definition a pre-protostellar object has to be gravitation- 
ally bound. In order to investigate the dynamical state of the 
core, we consider the virial theorem. The condition for a grav- 
itationally bound molecular clump or core is given by 



Emag ~t~ E po t 



2 ' (-Efci 



E K 



(3) 



where E mag denotes the net magnetic energy, E pot is the gravi- 
tational energy, Eki n the total kinetic energy and E ext accounts 
for external pressure. We assume that kinetic energy consists 



Mass 


50 ± 10 


M 


Dust Temperature 


12 ±2 


K 


Size 


0.55 x 0.29 


pc 2 


Visual extinction Ay 


30 


mag 


Column density n(H2) 


2.2 -10 22 


cm -2 


Density N(H 2 ) 


~ 2 -10 5 


cm" 3 


Gravitational energy 


6.4 -10 37 


J 



of the thermal energy of the molecules and the energy of their 
turbulent motion: 



Eh 



-NkT - 



-Ma 



turb 



(4) 



2 2 

We derive a\ urh = AV 2 /{Sln2) - kT/m, where AV is the 
observed line width of our ammonia data and m the mass of the 
ammonia molecule. The total line width of AV = 0.955 km s _1 
is an upper limit, since it was derived from a gaussian fit of the 
whole maingroup. On the basis on the (1,1) spectrum presented 
in Fig. H it is unclear if the two peak components correspond 
to the resolved hyperfine structure transition, which would in 
term lower AV, or if the double peak is simply due to noise in 
the spectrum. The error of the kinetic temperature Tki n =21.3 
± 5 K is rather high, due to the low SNR in the (2,2) spectrum. 
Usually the kinetic gas temperatures in cold and dense regions 
are lower than the obtained dust temperatures. We therefore 
consider Tkin as an upper limit. 

With the effective core radius of R = 0.2 pc we derive a 
gravitational potential of the homogenous ellipsoid (Liljestrom 
1991) 

E pot = 3GM 2 /5R (5) 

We find E pot = -6.4 • 10 37 J, E therm = 1.0 • 10 37 J and 
Eturb = 2.5 • 10 37 . Relative to \E pot \ the total kinetic energy 
sums up to 52 %. Neglecting external pressure and magnetic 
energy for which we have no data at present, we conclude that 
the cold core is roughly in virial equilibrium. If it will collapse 
or not, either as a whole or in fragments, will depend mainly 
on the local turbulent velocity field. Further molecular line data 
of sufficiently high spatial resolution are needed to assess the 
collapse issue. Tab. [7] summarizes the properties of the cold 
core region. 

4.2. The nature of the two cold submillimeter 
sources SMM1 and SMM3 

Phenomenologically, the youngest protostars (Class objects) 
are characterized by three attributes (Andre et al. 1993): 

1. Very little emission at A < 10/im, 

2. Spectral energy distribution similar to a blackbody of T ~ 
15-30 K and 

3. L submm /L bo i > 5 x 10~ 3 , where L submm is the luminos- 
ity measured at A > 350 /im. 

We suggest that SMM1 and SMM3 fullfil these criteria. SMM1 
and SMM3 are not detected in the mid-infrared, making them 
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Fig. 11. Spectral energy distribution of the compact submillimeter sources SMM1 and SMM3. The temperature is only indicated 
for = 1.5. See text for a discussion of the value of the /3-parameters. The arrow indicates a 60 /im upper limit based on IRAS 
photometry. 



candidates for the youngest protostellar objects in the vicinity 
oflSOSS J 20298+3559. 

The spectral energy distributions of SMM1 and SMM3 are 
presented in Fig. |TTJ Relying only on our submillimeter pho- 
tometry, the spectral energy distributions of both sources can 
be described by modified blackbodies with dust temperatures 
between 10 and 38 K. The large uncertainty of the temperature 
is mainly due to the unknown dust emissivity index (3 which 
covers the range between 1 .0 < (3 < 2.0 for protostellar objects 
and YSOs (Dent et al. 1998). To unambiguously determine the 
dust emissivity /?, measurements at the maximum or Wien part 
of the energy distribution are required. 

Next we try to provide 60 /^m upper limits for SMM1 
and SMM3, considering the flux of FIR1 minus that of main 
contributors. Due to confusion with large amounts of dis- 
tributed cold dust detected by the large-beam IRAS 100 /im 
and ISOPHOT 170 /im bands, the photometry at these wave- 
lengths can hardly be used for a reliable characterization of the 
two faint compact sources. The most stringent upper limit for 
SMM1 and SMM3 can be derived from the IRAS 60 /im band. 
The photometry at this wavelength is however affected by the 
emission of very small dust particles, transiently heated in the 
reflection nebula of the Herbig B2 star IRS 1 . The presence of 
such extended PAH and VSG emission is indicated on our mid- 
infrared images of IRS1. The IRAS measurements at 12/im, 
25/im and 60 /im are particularly sensitive to such an extended 
emission component. In order to determine the contribution 
of thermal radiation of transiently heated PAH and very small 
grains in these bands we have used the three-phase dust models 
by Desert et al. (1990). After subtraction of compact emission 
from IRS 1 according to our model (see next section), we de- 
termined the relative contributions of PAH and VSG emission 
from the residual IRAS 12/im and 25/im photometry and us- 
ing the model by Desert et al. (1990) for the radiation field of 
a stellar source with T e / / = 20 000K. We find a contribution 
of Fgo/OTi = 7 ± 3 Jy from the Herbig star IRS 1 and its reflec- 
tion to the total flux density of F e o^ m = 14 Jy. The residual flux 



density of Feofim = 7 Jy is considered as an upper limit for each 
of the two sources. 

Using the possible range 1.0 < (3 < 2.0 for the emissivity 
index, we derive total masses of 3.7 M < M < 18 M Q for 
SMM1 and 2 M Q < M < 8 M for SMM3. The corresponding 
bolometric luminosities cover the range 36 Lq > ^boi 
1.8 L for SMM1 and 34 L > Uoi > 3.3 L for SMM3 
respectively. If we favor an intermediate emissivity index 
of (3 ~ 1.5, which was the average value found in the large 
submillimeter survey of 73 young stellar objects by Dent et 
al. (1998), then we obtain for the two sources the following 
masses and luminosities: 
SMMl:M = 8.OM ,L 6oi =7L ,T d = 18 K 
SMM3: M = 3.3 M , L 6oi = 17 L , T d = 21 K 

We find L su b mm I L(, ; <; 1 % for each choice of 1 .0 < (3 < 
2.0, making SMM1 and SMM3 bona-fide candidates for Class 
objects. Further evidence for SMM1 and SMM3 being gen- 
uine protostars could be provided via detections of a molecular 
outflow and/or a cm-wavelength continuum source. 

4.3. Evidence for accretion onto the extremely 
young Herbig B2 star IRS1 

The Herbig B2e star IRS1 is the most luminous young stel- 
lar object embedded in the cold cloud complex. A very early 
evolutionary stage of IRS1 is indicated from its locus in the 
Hertzsprung-Russel diagram presented in Fig. [l~2|. According to 
the pre-main-sequence evolutionary tracks by Palla & Stahler 
(1993) the source is located close to the stellar birthline com- 
puted for M = 10~ 5 M yr _1 corresponding to a M star = 6.5 
M star with an age of less than 40000 yr immediately after 
becoming optically visible. 

The observation of inverse P Cygni profiles as shown by 
assymetric red and blue line components of the H I Balmer 



transitions presented in Fig. 13 provides evidence for ongoing 



mass infall in a magnetospheric accretion process (Edwards et. 
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Fig. 12. H-R diagram with stellar effective temperatures and 
luminosities, adopted from Palla & Stahler (1999). Stellar 
masses along the zero-age main sequence (ZAMS) are indi- 
cated in units of M Q . The Herbig B2 star IRS1 is located close 
to the stellar birthline, which was calculated for an accretion 
rate of M = 10 -5 M Q yr _1 . The pre-main-sequence tracks 
overlaid indicate a stellar age of less than 40000 yr. 



al 1994). From the redshifted absorption component in the Si II 
lines at 6347.1 A and 6371.4 A transitions shown in Fig. ^ we 
derive an infall speed of 1 10 ± 10 km s _1 . 

Near-infrared photometry and our sub-arcsecond resolution 
mid-infrared images revealed the presence of strongly con- 
densed (r ^ 250 AU) warm dust in the immediate surround- 
ing of IRS1. In order to model the optical and infrared spec- 
tral energy distribution of the heavily reddened source IRS1, 
a reliable determination of the extinction towards the source 
is essential. Based on the effective stellar temperature derived 
from our observations of photospheric He I lines, we apply a 
Kurucz model atmosphere with T e f f = 20 000 K for the stel- 
lar photosphere. We assume solar metallicity. Our dereddened 
BVRI photometry is best described for Ay = 6.1 ± 0.3 mag 
where the error is due to the uncertainty of T e ff ■ Assuming the 
distance d = 1.8 kpc of the cloud complex we derive a stellar 
radius of R* = 4 ± 0.5 Rq from the surface flux density of 
the Kurucz model. The integrated bolometric luminosity of the 
photosphere is L* = 2200 ± 600 L©. 

How can our observations of ongoing accretion and cir- 
cumstellar dust for this extremely young stellar object be com- 
bined in a consistent model? Fig. [l4] shows the dereddened 
spectral energy distribution of IRS 1 . The near-infrared emis- 
sion bump seen at about 3 /im wavelength is characteristic for 
Herbig AeBe stars and has been interpreted either by the pres- 
ence of a circumstellar disk with an inner hole (Hillenbrand 
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Fig. 13. Emission line profiles of the H I Balmer series. Ha, H(3 
and H7 are shown from top to bottom. The stellar continuum 
for a standard star of identical spectral type is shown for each 
transition. In addition to a strong central self absorption, the 
redshifted emission line component is suppressed for the higher 
order lines, indicating infall in an accretion process. 

et al. 1992), by spherical envelopes (Hartmann et al. 1993) or 
the combination of both (Miroshnichenko et al. 1999). While 
there is general consensus that the Herbig Ae stars are sur- 
rounded by circumstellar disks similar to those of T Tauri 
stars during most of their pre-main-sequence phase (Natta et 
al. 2001), only one out of seven Herbig Be stars was detected 
by mm-interferometry (Natta et al. 2000). Extended halos and 
envelopes in the near- (Leinert et al. 2001) and mid-infrared 
(Prusti et al. 1994) are more frequently observed around more 
massive Herbig Be stars. However, these observational differ- 
ences could be due to a faster evolution of massive stars and 
differences in the relative timescales of the PMS phases, which 
limit the detection of disks around early B stars to the very early 
stage of their evolution. 

We suggest that our spectroscopic signatures of ongoing 
accretion and the infrared excess can best be combined in 
terms of an accretion disk around IRS1. We have calculated 
the spectral energy distribution of an active viscous accre- 
tion disk including radiative heating by photospheric radia- 
tion from IRS 1 . The radial temperature dependence of the op- 
tically thick, geometrically thin disk is then assumed to be 
T(r) = [T QCC (r) 4 + T I . Qd (r) 4 ] 1 / 4 . We use a standard a-disk 
(Lynden-Bell & Pringle 1974). Its accretion disk temperature 
profile is given by 



Tacc(r) 



3GM*M 




g 7rcrr 3 



1/4 



(6) 



In the flat disk approximation its temperature profile due to ir- 
radiation only is 



T r ad{r) 



37rr 3 



-1 1/4 



(7) 
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Fig. 14. Dereddened optical and infrared spectral energy distri- 
bution of the Herbig star IRS 1 . We have modelled the strong 
infrared excess emission by an optically thick accretion disk. 
The star is described by a Kurucz model with a T e t /=20 000K 
photosphere (log g=4.0). Our observations obtained with the 
MPIA 0.7m telescope in the optical and with MAX at UKIRT 
in the thermal infrared are in good agreement with the model. 

Table 8. Properties of the Herbig B2 star IRS 1 and its accretion 
disk. 



Mass 


6.5 ± 1 


M 


Luminosity 


2200 ± 600 


L 


Effective temperature 


20000 ± 2500 


K 


Radius 


4.0 ± 0.5 


R© 


Visual extinction Av 


6.1 ±0.3 


mag 


Age 


< 40000 


yr 


Rotational velocity 


158 ± 25 


km/s 


Inner disk radius 


60.0 


R© 


Accretion rate 


5 x 10~ 5 


M yr 1 



Emerging spectral energy distributions have been com- 
puted for varying inner disk radii R^ and accretion rates M acc . 
As shown in Fig. [l4] the fit of the star+disk SED to the dered- 
dened photometry of IRS 1 is consistent with a dusty, optically 
thick accretion disk. Tab. |8| summarizes the model parameters 
of the system. 

Next we consider whether our data could also allow an en- 
velope: Our spectrophotometry of IRS 1 obtained in the narrow- 
band silicate filter at 9.7 /im shows a weak absorption fea- 
ture (Tg.j^m ~ 0.2 ± 0.1) which is consistent with the total 
line-of-sight extinction towards the source according to Whittet 
el al. (1987). Models of spherical envelopes which account 
for the observed spectral energy distributions of Herbig Be 
stars all predict strong silicate emission (Berilli et al. 1992, 
Miroshnichenko et al. 1997). In order to explain the suppres- 
sion of the silicate emission feature in a B0 star, the inner ra- 
dius of the envelope must be about 0.02 pc (Natta et al. 2000). 
This is inconsistent with our diffraction-limited mid-infrared 
images, which show that the bulk of flux from IRS 1 arises from 
a true point source, i.e. the radius of the mid-infrared emitting 
region is < 400 AU. 



4.4. Star formation efficiency 

We estimate the star formation efficiency SFE = M star s I M tot 
in the central region ISOSS J 20298+3559-FIR1 . M tot includes 
the masses of the diffuse component, the dense core and all 
young stellar objects. We have determined the masses for the 
near-infrared objects IRS2..5 from the absolute magnitude- 
mass relation by Carpenter et al. (1993). 

Concerning the masses of the Class sources SMM1 and 
SMM3 we assume that they loose 50 % of their envelopes dur- 
ing their further evolution (Bachiller et. al. 1996). The mass of 
all young stellar objects is then M sta rs= 20 M . Adding the 
mass of gas and dust in FIR1 to the stellar masses we obtain 
Mtot = 140 M Q . This results in a star formation efficiency of 
about 14 %. 

To compare: The SFE is about 2-9 % for the Taurus com- 
plex (Cohen & Kuhi 1979) and 1-9 % in Orion A+B (Carpenter 
2000). ISOSS J 20298+3559-FIR1 clearly has at least the effi- 
ciency of these low- and high-mass star forming regions. 

5. Conclusions 

We have identified the young star forming region ISOSS J 
20298+3559 performing a cross-correlation of cold compact 
far-infrared sources from the ISOPHOT 170 /im Serendipity 
Survey database with the 2MASS, MSX and IRAS surveys. 
Multi-wavelength follow-up observations of this region yield: 

1. The star forming region is associated with a complex of 
four optical dark clouds C1..C4 which have a total mass of 
760 Mq. 

2. We derived a distance of 1800 ± 300 pc based on op- 
tical extinction data. This associates the region with the 
CygnusX Giant Molecular Cloud in agreement with our 
molecular line kinematics. 

3. The cold ISOSS source FIR1 corresponds to the dense in- 
ner region of the central dark cloud C 1 and contains a total 
mass of 120 M Q gas and dust with an average temperature 
ofl6K. 

4. We have identified two candidate Class objects SMM1 
and SMM3. The sources have masses of 8 and 3.5 M Q 
which makes them precursors of intermediate mass stars. 

5. The externally heated cloud core of CI has a total mass 
of 50 M Q and a central dust temperature as low as 11K. 
Ammonia in the NH3(1-1) transition has been detected. 
The object is gravitationally bound as derived from our am- 
monia molecular line observations, which makes it a can- 
didate massive pre-protostellar core. 

6. The most luminous object in the vicinity is the Herbig B2 
star IRS1. The object has a mass of 6.5 M and a bolomet- 
ric luminosity of 2200 L Q . Inverse P Cygni profiles in the 
higher HI Balmer series and Sill lines indicate ongoing ac- 
cretion. The spectral energy is well described by a Kurucz 
model for the stellar photosphere and a viscous reprocess- 
ing accretion disk. The stellar age inferred from pre-main- 
sequence evolutionary tracks is less than 40000 yr. 

7. Several embedded near-infrared sources have been identi- 
fied. One of them (IRS7) is surrounded by a mid-infrared 
reflection nebula. 
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8. The star formation efficiency in the dense and cold region 
FIR1 is about 14 %. 
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